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We propose a general formalism to compute exact correlation functions for Cardy's boundary 
states. Using the free-field construction of boundary states and applying the Coulomb-gas technique, 
it is shown that charge-neutrality conditions pick up particular linear combinations of conformal 
blocks. As an example we study the critical Ising model with free and fixed boundary conditions, 
and demonstrate that conventional results are reproduced. This formalism thus directly associates 
algebraically constructed boundary states with correlation functions which are in principle observable 
or numerically calculable. 
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I. INTRODUCTION 

Modular invariance of partition functions plays ex- 
tremely important roles in two-dimensional conformal 
field theory (CFT). The ADE classification of modu- 
lar invariants by Cappelli, Itzykson and Zubcr [1,2] (sec 
also [3-5]) is obtained by considering CFTs on the torus. 
The classified modular invariants correspond to particu- 
lar sets of operators, which are supposed to model critical 
systems in certain universality classes. Similar consider- 
ation also applies to CFTs defined on a manifold with 
boundary. For a CFT on the cylinder, the constraints 
from modular invariance lead to a classification of bound- 
ary states. This method was invented by Cardy [6] in the 
eighties and developed by many in the nineties [7-12]. 
Such a classification of boundary states with consistent 
modular properties has recently attracted much attention 
along with the development of D-brane / string theory 
and various application of boundary CFT to statistical 
physics. 

Boundary states with consistent modular properties, 
or consistent boundary states, are normally believed to 
represent boundary conditions which may be physically 
imposed on D-branes or borders of statistical systems. 
In order to understand the behaviour of observables in 
the presence of such a boundary, we need to find cor- 
relation functions for consistent boundary states which 
are defined through Cardy's classification. In principle 
this can be accomplished by operator product expan- 
sions (OPEs), that is, by finding boundary operators for 
a given boundary state, solving the constraints satisfied 
by coupling constants, and then obtaining the correla- 
tion functions by OPEs using the boundary operators. 
The correlation functions obtained in this way are per- 
turbative, that is, in the form of series expansion. For 
practical use, we often need to know exact correlation 
functions and in such a case we normally solve differ- 



ential equations to find conformal blocks and fix their 
coefficients by physical considerations [13]. In this differ- 
ential equation method, however, the relation to Cardy's 
classification of boundary states is not quite evident. 

In this paper, we present an alternative formalism 
for finding boundary correlation functions* directly us- 
ing boundary states obtained by Cardy's method. It is 
well known in string theory that correlation functions 
are simply given by inserting vertex operators within 
amplitudes (with or without boundaries). This picture 
is generalised to non-bosonic (c ^ 1) CFTs using the 
Coulomb-gas formalism on arbitrary Riemann surfaces 
(without boundaries) [14-17]. The purpose of this paper 
is to present a method to calculate boundary correlation 
functions based on the Coulomb-gas picture, by using the 
free-field representation of boundary states developed in 
[18]. Driven by a similar motivation, Coulomb-gas sys- 
tem on the half plane is discussed in [19], where the Ising 
model conformal blocks are reproduced using the contour 
integration technique of [20] and insertion of boundary 
operators. The key ingredient of the formalism proposed 
in the present paper is boundary states having bound- 
ary charges, which account both for the contour integral 
expression of conformal blocks and for their coefficients. 
The role of such boundary states has not been fully inves- 
tigated so far in this context, and our formalism allows 
a systematic study on the relation between correlation 
functions and algebraically defined boundary states. In 
the following we shall mainly consider Coulomb-gas sys- 
tems on the unit disk, where Felder's charged bosonic 
Fock space (CBFS) construction [14] is readily used. Al- 
though our method itself is quite general, we shall focus 



*We use this term for correlation functions of bulk operators 
1 the presence of boundary. 
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on presenting the ideas in simplest cases and show that 
it reproduces known results obtained by the conventional 
approach. 

We organise the rest of this paper as follows. In the 
next section we fix our notations and review the free-field 
construction of boundary states [18] . We describe in Sec. 
Ill our method to compute correlation functions on the 
disk and on the half plane. In Sec. IV we illustrate the 
method using the Ising model and show that it repro- 
duces the results of [13]. Finally in Sec. V we summarise 
and conclude. 

II. BOUNDARY STATES IN COULOMB-GAS 
FORMALISM 

Let us start, for the sake of self-containedness, by 
summarising the Coulomb-gas construction of boundary 
states [18]. The idea is to define coherent states in the 
charged bosonic Fock spaces (CBFSs) and find conditions 
for their diffeomorphism invariance and modular consis- 
tency. 

A. Coulomb-gas and charged bosonic Fock space 

In the Coulomb-gas formalism [21], the Virasoro min- 
imal models arc realised by the action, 

j d'x^id^^d^'^ + 2^/2a(ii^R), (1) 

where the scalar field <^{x) is assumed to decouple into 
two chiral sectors, ^{z,z) = (p{z) + <p{z). The energy- 
momentum tensor is obtained by the variation of the ac- 
tion (1) as 

T{z) = = - ^ : dipdip : +iV2aod'^ip. (2) 

The central charge of this system is 

c = 1 - 24al. (3) 

The vertex operators, 

Va{z) =: e'^°^(^) :, (4) 

are chiral fields of conformal dimensions ha = — 2aoa. 
In particular, the primary fields (j)r.s (0 < r < p', < s < 
p) of a minimal model are realised by the vertex operators 
Var,s (z) with 

ar,s = ^{^-r)a+ + ^{l- s)a_, (5) 

where a-|- — \/p/p', a- = —\/p'/p, and p and p' (we 
assume p> p') are the two coprime integers characteris- 
ing the minimal model. The conformal dimensions of the 



operators are 

hr,s ^ '^{ra+ + sa-)'^ - al, (6) 

which agree with the Kac formula. 

The holomorphic chiral boson field is expanded in 
modes as 

(p(z) = (fo - iaolnz + —z'", (7) 

and likewise for the antiholomorphic counterpart, 

(p{z) = (po — iao \nz + i —z~". (8) 

As we try to consider boundary CFT in the Coulomb-gas 

picture, a subtlety arises in the treatment of zero-mode, 
since the zero-mode of ^{z, z) does not naturally decou- 
ple into the holomorphic and antiholomorphic sectors. In 
our formalism, we shall simply split it into two identical 
and independent copies. In this sense, the boundary the- 
ory we shall consider is not exactly a non-chiral theory 
on a manifold with boundary but rather a chiral the- 
ory on its Schottky double^. We thus have two copies of 
Heisenberg operators, satisfying the algebra 

[Om, On] = rndm+nfi-, [fO, «o] = 

[am-, an] = mSm+n,o, [fio, Oo] = i, (9) 
with no interaction, 

[am, an] = 0, [(fio, ao] = [<Po, ao] = [vo, <fio] = 0- (10) 

In terms of the Heisenberg operators, the Virasoro oper- 
ators are written as 

i„^o = ^X^fln-feafe - \/2ao(n+ l)an, (H) 

fcez 

Lo = ^2 a-kak + ~ V^aoao, (12) 
fe>i 

and likewise for the antiholomorphic counterparts. 

The Hilbert space of CFT is reahsed in CBFSs with 
ERST projection [14]. For a chiral theory the CBFS 
Fa-ao is defined as a space obtained by operating with 
the creation operators a„<o on the chiral highest weight 
state ]a;ao) (see [14,18,23]). Since the boundary inter- 
twines the two chiral sectors, we need to construct non- 
chiral Fock spaces J^-'q,,-,;,!,, in order to describe boundary 
states. We denote the Mobius invariant non-chiral vacua 



^The Schottky double is a Riemann surface obtained by dou- 
bling the manifold except for boundaries. See, e.g., [22]. 
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with background charge ao as^ (0, 0;q;o| and |0, 0;q:o), 
and define highest weight vectors {a, d; ao\ and \a, d; ao) 
as 



(a,a;ao| = (0, 0;ao|e" 



i\/2a(po p—iV2a(po 



\a, d; ao) = e'V^«*'°e'^"'^° |0, 0; 



(13) 
(14) 



The state (a,a;ao| has holomorphic and antiholomor- 
phic charges —a and —a, respectively. Likewise, 
|a, a; ao) has holomorphic and antiholomorphic charges 
a and d. Using the Heisenberg algebra (9) it is easy to 
verify that these states satisfy 



(a, a; ao|ao = (a, a; ao |V2a, 
(a, a; ao|ao — (cK) ck; ao|\/2a, 
ao\a,d;ao) = \/2a|a, a; ao), 
ao\oe, a; ao) = \/2a|a, a; ao). 



(15) 
(16) 
(17) 
(18) 



These highest weight vectors are eigen states of the Vi- 
rasoro zero-modes: 

Lo|a, a; ao) = (a^ - 2aao)|a, a; ao), (19) 
Zo|a, a; ao) = (a^ — 2aao)|a, a; ao). (20) 

The states \a, d; ao) are annihilated by a„>o, an>0) Ln>Oi 
Ln>a, and (a, a; ao| are annihilated by a„<o, a„<o, Ln<Oi 
Ln<o- The non-chiral CBFS J-'a,a;ao is built on the high- 
est weight vector |a, a;ao) by operating with a„<o and 
a„<o. The dual space ^a.a-.ao defined similarly, by 
acting with a„>o and a„>o on (a, a;ao|. The non-chiral 
CBFSs thus defined are essentially the direct products of 

Chiral CBFSs, J^a,a;ao = Fa-ao ® Fa-ao- 

The inner products of highest weight vectors are sub- 
ject to charge neutrality, i.e. they are non-vanishing only 
if the net charges in the two sectors are both zero. The 
normalisation of the highest weight vectors must be in 
accordance with this condition and thus we have 

(a, a; ao|/3, ^; ao) = K^a./jJ^ j. (21) 

In particular, the vacua arc normalised as 

(O,0;ao|0,0;ao) = AC. (22) 

In unitary theories the constant k is usually positive and 
we normalise it to unity. As the Coulomb-gas system may 



■''Here wc give the same background charge ao to both sec- 
tors. Even if we relax this condition and start by allocating 
different background charges ao and ao to holomorphic and 
antiholomorphic sectors, respectively, the condition (26) re- 
stricts either ao = ±ao. For ao = — ao we have Q, = —1 and 
a — a, — 2ao = instead of (29) and (30). This merely flips 
the sign of all antiholomorphic charges and thus does not give 
any new results. 



well include non-unitary theories, k can be negative. In 
that case we choose k — 1. Thus we have 



u{a, a; ao\(3, (3; ao)u = Sa,p5^ p, 
N{a,d;ao\P,(3;ao)N = ^Sa,pS^p. 



(23) 
(24) 



The subscripts U and N stand for unitary and non- 
unitary sectors, respectively. Those two sectors have no 
intersection. 



B. Diffeomorphism invariant boundary states 

Boundary states appearing in CFT are diffeomorphism 
invariant in the following sense [6,24]. Let us consider a 
CFT on the upper half plane Im^ > 0, where C is a com- 
plex coordinate, = x + iy, x, y € The boundary is 
y = 0, or C = C*. Since the antiholomorphic coordinate 
dependence ( on the upper half plane may be mapped 
into the holomorphic dependence C* on the lower half 
plane [13], we often identify ( with (* . Now, once the 
boundary is fixed, the conformal symmetry of the the- 
ory should be restricted so that the boundary is kept 
fixed. In other words, as the conformal transformation is 
generated by the energy-momentum tensor, the energy- 
momentum flow across the boundary must vanish. 



[T{0-f{0] . = 0. 



(25) 



This may be translated into a condition on boundary 
states by mapping a scmiannular domain on the upper 
half C-plane into a fuU-annulus on the z-plane by z = 
exp(-27ri^/L), ^ = (T/7r)lnC. The boundary of the C" 
plane is mapped to the two concentric circles \z\ = 1, 
exp(27rr/L) bordering the annulus on the z-planc. Since 
the z-plane allows radial quantisation, (25) is written as 
the condition on the boimdary states \B) (on |2:| = 1), 



{Lk-L_k)\B)=0. 



(26) 



This condition, often called the Ishibashi condition, must 
then be satisfied by any boundary state in CFT. 

We may follow the standard construction of bound- 
ary states in open string theory [25-27] and find bound- 
ary states on CBFS by starting from the coherent state 
ansatz. 



K(-Ba,a;ao;n| = kI^, a; ao| exp ( 



fc>0 



afeOfc 



(27) 



\Ba,a;ao;n)K = ]J (=xp f-^a_fca_fc j |a,a;ao)K. 
k>0 ^ ^ 

(28) 

The subscript k is either U or N, specifying the normali- 
sation of the vacuum. As we have expressions of the Vira- 
soro operators (11), (12) written in terms of the Heisen- 
berg operators, one can see how the Virasoro modes L„ 
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and Z_„ operate on the coherent states \Ba.a:ao:ii) k by 
an explicit computation. It is shown [18] that the condi- 
tion (26) is satisfied if 



and 



= 1, 



a + a — 2ao = 0. 



(29) 



(30) 



Similarly, we see that K{Ba,a;ao;n\{Ln - L^n) = as 
long as (29) and (30). In the following we shall only con- 
sider such manifestly diffeomorphism invariant boundary 
states satisfying (29) and (30), and for simplicity we de- 
note, 



K{B{a)\ — K(-Ba,2ao-a;ao;f^=l|! 
|-B(Q;))k = \Ba,2cxo-a;ao;n=l). 



(31) 
(32) 



We note that the sum of the (holomorphic + antiholo- 
morphic) boundary charges agrees with the topological 
background charge on the Schottky double. Due to the 
condition (30), an inner boundary (on the z-planc) c;on- 
tributes 2ao to the sum of the charges, and an outer 
boundary contributes —2ao. When we consider an an- 
nulus whoso Schottky double is a torus, the sum of the 
boundary charges is zero (2ao'~2ao = 0), which coincides 
with the background charge of the torus expected from 
the Gauss-Bonnet theorem (the Euler number of a torus 
vanishes). For a disk, there is only an outer boundary 
which gives a charge —2ao. This agrees with the back- 
ground charge of a sphere, which is the Schottky double 
of the disk. Thus, the geometry of bulk can be assumed 
to be flat everywhere, since the curvature of the manifold 
is concentrated on boundary. 



C. Ishibashi states in free-field representation 

As the basis of boundary states in CFT is nor- 
mally spanned by Ishibashi states, wc need to construct 
Ishibashi states in terms of our Fock space representation 
in order to translate existing results of boundary CFT 
into Coulomb-gas language. At least for the diagonal 
minimal models, it is shown [18] that the Ishibashi states 
are expressible as linear combinations of [^(a))^, as far 
as partition functions on the cylinder arc concerned. 

Ishibashi states are defined for chiral representations of 
CFT and diagonalise the cylinder amplitudes (overlaps) 
to give characters. 



~l/2NLo4-io-c/12 



(33) 



Here, we are considering a cylinder of length T and cir- 
cumference L, or equivalently, an annulus on the 2:-plane 
with 1 < l-zj < exp(27rr/L). As the Hamiltonian is writ- 
ten as H = {2iT / L){Lq + Lq — c/12), the left hand side of 



(33) is ((ije-^^jj)) with q = er^'^'^/^. The characters of 
the minimal models are given by Rocha-Caridi [28], 



X(r,s)(a) 



Tr(.,«)<?^°-^/'' 

e 



pr—p's^pp' 



(t) ^pr-{-p' s,pp' ('^) 



(34) 



where Q\,fj,{T) and r]{T) are the Jacobi theta function 
and the Dedekind eta function, defined as Qx,h{t) = 
^^^^g(2^fe+A)V4M and r?(T) ^ q'/'''Un>li^ " with 
q = e^'^*'^. Thus, the cylinder amplitudes (33) are power 
series in q, divided by ri{f). 

As our boundary states \B{a))K are defined in a Fock 
space representation, we may explicitly compute the 
cylinder amplitudes between such states. They are [18], 

.{B{a)\e-^"\Bm, 
= .{B{a)\iq'/')'^°+^°-'^/'^\Bm, 



^(a-ao) 



-kS, 



(35) 



The amplitudes between unitary and non-unitary sectors 
(e.g. u {B {a)\e^'^ ^ \ B ( P)) n) vanish because these sectors 
do not intersect. Note that ^("""o) /r]{T) is the character 
Xa;ao {q) of the chiral CBFS Fa-a^ . In this sense, the state 
\B{a))f^ may be regarded as the Ishibashi state of Fa-ao- 
We may now compare the expressions (33) and (35) to 
find a possible free-field representation of the Ishibashi 
states of minimal models. Defining 



and 



with 



|(r, s))) = \ar,s)u + \ar,-s)N, 



u{ar,s\ = X] u{B{ar,s + k^/pp')\, 

N{ar-s\ = ^N{B{ar-s + k^/pp')\, 

feez 

\ar,s)u = X] \^i'^r,s + k^/pp'))u, 

fcez 

\ar-s)N = ^ \B{ar,-s + ky/p^))N, 



(36) 

(37) 

(38) 
(39) 
(40) 
(41) 



feez 



where a^.s are given by (5), one can show that the states 
|(r, s))) diagonalise the overlaps and give minimal model 
characters, 

(((r,s)|(«~l/2)^«+^c-/12|(/,5'))) = 5rr'5ss'X(.r,s){q)- 

(42) 

The appearance of the non-unitary sector even in uni- 
tary CFTs may seem odd, but this is necessary to 
describe cylinder diagrams where otherwise unphysical 
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states would propagate. On the disk the non-unitary 
sector decouples and does not contribute to correlation 
functions (Sec. III). The states (36), (37) are a good 
candidate for the minimal model Ishibashi states as far 
as the modular properties are concerned. We, however, 
immediately notice that such 'Ishibashi' states are not 
unique since \(r,s)) and \{p' — r,p— s)) give rise to the 
same character but are perpendicular to each other. In 
order to have a unique Ishibashi state for each primary 
field (r, s) ^ (p' — r,p — s) of minimal models, we define 
the symmetrised states, 

((</'r,s| = {{(t>p'-r,p-s\ 

= ±{{{{r,s)\ + {{{p'-r,p-s)\), (43) 

|0r,s)) = |0p'-r,p-s)) 

= l={\{r,s))) + \{p'-r,p-s)))), (44) 

and shall regard them as our free-field realisation of the 
Ishibashi states. Although this 'symmetrisation' was not 
considered in [18], such a prescription to ensure the equiv- 
alence of (r, s) ~ {p' — r,p— s) is necessary. 

D. Cardy's consistent boundary states 

Physical boundary states in CFT are not only diffeo- 
morphism invariant, but must also satisfy an extra con- 
straint called Cardy's consistency condition. We consider 
a cylinder of length T and circumference L as before 
and assume that boundary conditions a and /? are im- 
posed on the two boundaries. Then, depending on how 
we define the direction of time, the partition function 
on this cylinder may be written in two different ways. 
We may first regard the cylinder as an open string prop- 
agating in the periodic direction of time, with bound- 
ary conditions a and /3 imposed at the two ends of the 
string. The partition function is then a sum of the chi- 
ral characters, Z^^{q) = T,j'"'^^^Xj{'l), where XjiQ) is 
the character for the chiral representation j and rv'^ - is 

afj 

a non-negative integer representing the multiplicity of 
the representations. We have defined q = e~'^^^^. We 
may also see the cylinder as a closed string propagat- 
ing from one boundary (with boundary condition /3) to 
the other (with d). Then the partition function is sim- 
ply the cylinder amplitude between the two boundaries, 
{a\e~'^^\/3). Due to the equivalence of the two pictures, 
we have T,jn{^Xjiq) = {d\{q^/'^)^°+^'>-''/'^^0}. This is 
the consistency condition which needs to be satisfied by 
the boundary states {d\ and \l3). 

If we have an appropriate basis of the boundary states, 
the right hand side of the consistency equation may be 



expanded using the basis states {{a\}, {\b)} as 

J2<^XM = J2i'^\a){a\{q'/Y''^^''-''^''\b)m. 

j a,b 

(45) 

Solving this equation, the consistent boundary states 
{la)} are found as linear combinations of the basis states. 
It is convenient to use the Ishibashi states for such a 
basis. For diagonal minimal models, using the prop- 
erty (33) of the Ishibashi states and the modular trans- 
formation of the characters Xii^) = J2j ^ijXjiQ) under 
T ^ f = — 1/r, we have, by equating the coefficients of 
the characters, 

J2ni^S,j = {a\j)){m. (46) 

i 

Assuming the existence of a state |0) such that rig. = 
nig = (5g, (46) was solved by Cardy [6] as. 

Now, as the minimal model Ishibashi states have been 
written in the free- field representation (43) (44), Cardy's 
consistent boundary states (47) can also be expressed 
using our coherent boundary states by substituting (44) 
into (47). 

For the convenience of later discussions, let us spell 
out these in the specific example of the Ising model. The 
Ising model is the simplest non-trivial minimal model 
having the central charge c = 1/2 and is characterised 
by the two coprime integers p = A and p' = Z. There are 
three operators, the identity /, the energy e and the spin 
cr, having the conformal dimensions 0, 1/2, and 1/16, 
respectively, and are identified in the Kac table as / = 
•^1,1 = 02,3, e = (/)2,i = 01,3 and a = (/)i,2 = 02,2- Using 
the modular transformation matrices for the characters, 
Cardy's boundary states are written as [6] 

|7) = |0) = 2-i/2|7)) + 2-i/2|e)) + 2-^/^\cj)), (48) 
|e) = 2-^'^\I)) + 2-V2|e)) - 2-^'^)), (49) 
|a) = |/))-|e)). (50) 

It is argued that the first two states correspond to the 
fixed (up and down) boundary conditions since they differ 
only by the sign of |cr)) which is associated with the spin 
operator. The last state (50) is then identified as the free 
boundary state. Using (37) and (44), these states are 
written in our free-field representation as 

\i) = 2"^(|ai,i)c/ + |ai,-i)Af + \a2.?i)u + \a2-z)N 
+ |a2,i)c/ + |a2,-i)Af + \ai,i)u + \ai,~z)N) 

-h2-3/4(|ai,2)c/ + |ai,-2)Ar + \a2,2)u + |a2,-2)Ar), 

(51) 
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|e) = 2"^(|ai.i)f/ + |ffli.-i).Y + |a2,3)c/ + 1^2, -3) 
+ |a2,l)(7 + |a2,-i)Ar + \ai,3)u + |ai,-3)Af) 

-2-3/4(|ai,2)[/ + |ai,_2)jV + |02,2)[/ + |02,-2>w), 

(52) 

= 2~^/^(|ai,i)i7 + iai ,-i)jv + \a2,3)u + ia2 ,-3)jv 

— \(l2,l)u — |a2,-l)jV — \0'l,3)u — |«l,-3)jv), 

(53) 

where the states on the right hand sides are defined by 
(40), (41). They are superpositions of countably many 
coherent states with diflferent boundary charges. 

III. BOUNDARY CORRELATION FUNCTIONS 

Now let us discuss how to compute boundary correla- 
tion functions in our Coulomb-gas picture. After giv- 
ing the general formalism, we shall focus on the one 
point function of (j)r,s{z,z) and the two point function 
of 4'i^2{z, z), and derive their explicit expressions on the 
unit disk. Once correlators on the disk are obtained, it 
is straightforward to map them on the half plane. At the 
end of this section, we compare our free-field approach 
and OPE computation of boundary correlation functions. 

A. Screened vertex operators and BRST states 

In the CBFS language, correlation functions on the 

full plane are described as follows [14]. We define chiral 
screened vertex operators V^g'"{z) as 

/m n 
n du, H dv,VrAz)V+{ui) ■ ■ ■ V+{Um) 
i=\ j = l 

xV-{Vi)---V-{Vn), (54) 

where for conciseness we have denoted Va^ ^ (z) as Vr.s{z) 
and Vaj.{z) as V±{z), and the integration contours are 
those of Felder's, all going through z and encircling the 
origin (Fig.l of [14]). Such an operator is a primary field 
of conformal dimension hr,s- We also denote the chiral 
CBFS i^c«,,,;ao as Fr,s- The operator K^'^C-^) defines a 
map from one Fock space to another, 

(55) 

The p-point correlator, 

(0; aoKZr {z,) ■ ■ ■ KT^;"''(^P)|0; «o), (56) 
is then seen as a sequence of mappings, 

Fl,l ^ Frp-2mp,Sp-2np ^ ' ' ' , (57) 

and the final state KT,y,"H^i) • • • Vr'^^s'"" {zp)\0; Qq) must 
belong to Fi^i in order to have a non-trivial inner product 



with (0;q:o| <= F* (the dual module of i^i,i). The same 
p-point correlator may be expressed differently, as 

ao|KT,'y"n^i) • • • VripHzp)\0; Qq), (58) 

but this is in fact proportional to (56). A key object in 
this formalism is the BRST operator, 

^27Tia^r f ^ 

Qr = —^2 f'[[dUiV+{u^)---V+{Ur), (59) 

which maps JV,s to F-r,s- The BRST operator is nilpo- 
tent, QrQp'-r = 0, and physical states are realised as the 
cohomology space, KeiQr/ImQp'-, - An important point 
which is evident in this picture is that all intermediate 
states appearing in the correlator are BRST states since 
the vacuum |0; ao) is a BRST state and the screened ver- 
tex operators map BRST states to BRST states [14]. 

We shall combine the above machinery and the bound- 
ary states of the previous section to compute correlation 
functions on the unit disk. The in-state of the correla- 
tors is the non-chiral vacuum (in the unitary sector) at 
the origin, [0) = |0,0;ao)c/, whic:h is a BRST state. For 
the out-state, we choose a boundary state u{B{a)\ of 
(31) with a fixed boundary charge a. We only consider 
the unitary sector since the non-unitary sector does not 
give non-trivial inner products with the in-state vacuum. 
The correlators are then obtained by inserting non-chiral 
screened vertex operators, 

^&"£f ^(^^'^) = (^ot^r (^), (60) 

between u{B{a)\ and |0,0;ao)c/- As we focus on di- 
agonal theories, the antiholomorphic vertex operators 
have the same conformal dimensions as the holomorphic 
counterparts, i.e., either (fj,Si) = {ri,Si) or {fi,Si) = 
(p' — ri,p— Si). Note that, in this construction, all the in- 
termediate states are manifestly BRST invariant because 
no spurious states arise. Actual boundary p-point cor- 
relation functions for physical boundary conditions are 
obtained by summing the fixed boundary-charge correla- 
tors, 

t/(B(a)| nK^';"H^i)V^'J:r^(^i)|0,0;ao)c/, (61) 

over the boundary charges according to the linear com- 
binations (such as (51) - (53)) given by Cardy's method. 
Note that (61) is non- vanishing only for certain config- 
urations of screening charges. The net holomorphic and 
antiholomorphic charges are respectively, 

— a + ^ ar,,si + ^ + ^ '^^Q!-, (62) 

i i i 

and 

a - 2ao + ^Q!f^,si + ^'7iia+ -h ^n^a-, (63) 

i i i 
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which must vanish independently. These charge neutral- 
ity conditions associate the allowed values of a with the 
numbers of holomorphic and antiholomorphic screening 
operators. 

In the computation of the expression (61), 

Vr,siz)Vr,s{z), Vr^siz)Vp' -r,p-s{z) , Vp' ^r.p~siz)Vr^siz) 

and Vp'-r,p-s{z)Vp'-r,p-siz) all correspond to a non- 
chiral field (f>r,s{z,z) and one may use whichever com- 
binations one wishes to use. This is ensured by the fact 
that (61) is essentially a chiral 2p-point function where 
the equivalence of (r, s) <-^- {p' — r,p — s) (after the trun- 
cation of unphysical states) is guaranteed [14,21]. In par- 
ticular, we are allowed to use V^,s(z) as the holomorphic 
and Vpi-r,p-siz) as the antiholomorphic ('mirror image') 
part of a single non-chiral operator (the analogy of a 
mirror and a boundary of CFT is based on the confor- 
mal Ward identity [13] which does not distinguish V5.,s(^) 
from Vp' -r.p-s{z))- Due to this equivalence, apparently 
different choices of vertex operators should all lead to a 
same conformal block function. In practice, likewise to 
the case of the Coulomb-gas computation without bound- 
ary, we shall choose such vertex operators that the num- 
ber of screening operators is minimised. In the following 
subsections we shall evaluate the expression (61) for par- 
ticular cases of one and two point functions. 

B. Boundary one point functions 

For evaluation of the boundary one point correlator, 

u{B{a)\V;:i'-{z)V,YimO;ao)u, (64) 

it is convenient to choose (f , s) — {p' — Vjp— s) (the other 
choice (f, s) — (r, s) should lead to the same result but 
involves complicated integral expressions). According to 
(55), the holomorphic part of the CBFS is mapped as 

Fis Fr~2m,s-2n, (65) 

and the antiholomorphic part is mapped as 

Fis — > Fp'-r-2m,p-s-2n- (66) 

Since u{B{a)\ e F*.^^ (g) F^^^_^.^^^, the correlator 
is non-vanishing only when Fa-ao = -Fr-2m,s-2ri and 

F2aQ~a-,aQ — -^p'— r— 27n,p— s — 2ni that is, 

a = i(l - r + 2m)a+ + i(l - s + 2n)a_, (67) 

and 

2ao~a^ l-{l-p'+r+2m)a+ + l-{l-p+s+2n)a^. (68) 



Summing (67), (68) and using a+ = \/p/p', a- = 
— y/p'/p, we have (to + m)a_|_ + {n + n)a- — 0, or 

TO-|-TO = 0, n + fi — 0, (69) 



implying no screening charges. Then from (67) we have 
a = ar^s, and the correlator (64) is evaluated as 

u{Biar,.)\V,y{z)V,'f{z)\0,0;ao)u - {I ~ zz)-'\ (70) 

where h — ar,s(Q^r,s ~ 2ao). As we shall see later, this is 
proportional to the two point correlator on the full plane. 



(a) Conformal block /j. 




(b) Conformal block In. 




FIG. 1. The two conformal blocks Ji and In in the bound- 
ary two point function of (f)i^2 on the disk. The sphere repre- 
sents the double of the disk, and the upper and lower hemi- 
spheres stand for the holomorphic and antiholomorphic sec- 
tors, which are glued on the boundary (the equator). The 
lower hemisphere coordinates z* are obtained from Zi by the 
doubling {zi, Zi) — > [zi, z*). In the case of Ii where the screen- 
ing operator Q- lie in the holomorphic sector, the integration 
contour can be deformed into the Pochhammer type around 
zi and Z2. The integral is then proportional to the one from 
Z2 to zi (a). Similarly, the screening contour of In is in the 
antiholomorphic sector and the integral is proportional to the 
one from to Zi (b). 
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C. Boundary two point functions of 01,2 

As a less trivial case, we consider the two point 

correlator of the primary field (t>i_2- For the conve- 
nience of calculation we define one of the operators as 
01,2(2:1,21) = Vi, 2(2:1) 14"-i,p-2(^i) and the other as 
f/*!, 2(^2, ^2) = Vi,2(2;2)Vi,2(^2)- The expression (61) then 
becomes 



u{B{a)\v,y'-^{z,)v;rx;_,{z,) 

xyi7'"^(^2)n7"*^(^2)|0,0;ao)c/. 



(71) 



In the holomorphic and antiholomorphic sectors, the 
CBFSs are mapped as 

-F'l.l — * -F'l-2m2,2-2n2 ~^ -Fl-2mi-2m2,3-2ni-2n2 i (^2) 
-Fl.l — * -Fl-2m2,2-2n2 ~^ -Pp'-2mi-2m2-l,p-2fii-2ra2-l) 

(73) 

and hence, in order that the correlator be non- vanishing 

we must have 

a = {mi + m2)a+ + (ni + n2 - l)a-, (74) 
2ao - a = (mi +m2 + l)a+ + (ni + n2 + l)a-. 

(75) 

Adding the above two expressions we have (m + m)a+ + 
{n+n—l)a- = 0, where m = mi+m^ and n = ni+n2 are 
the numbers of positive (0;+) and negative (a-) screen- 
ing operators in the holomorphic sector, and similarly 
m = mi + m2 and n = ni+n2 for the antiholomorphic 
counterparts. This charge neutrality condition implies 



m-|-m = 0, n + n = l. 

Then we have two possibilities: 

(I) m = ffi = n = 0, n= 1, 
(II) m = fh = n = 0, n= 1. 



(76) 



(77) 
(78) 



From (74), we have the boundary charge a = ai,i = 
for (I) and a = ai,3 = — a_ for (II). 

For the charge configuration (I), we have one screening 
operator 



Q- = j>dvV-{v), 



(79) 



in the holomorphic sector and thus the correlator (71) 
takes the form. 



,l)\fdvVi,; 

xVi,2{z2)Vi,2{z2)\0,0;ao)u, (80) 
which, by an explicit calculation, reduces to 



^ dv{l - ZiZiy{l - ZiZ2)\l - VZiYil - VZ2) 



x{l-Z2ZiY{l-Z2Z2)\zi- 
X{Z1 - Z2f{v - Z2f{zi - Z2) 



(81) 



with a = 2ai,2(2ao — 0:1,2), b = 2ai 2, c = 2a_(2ao — 
0:1.2)5 ^iiid d = 2q!_q:i.2- As this expression is analytic, 
we may deform the integration contour as long as it is 
closed and non-contractable. In this case the screening 
operator must lie entirely on the holomorphic part and 
Felder's contour can be deformed into the Pochhammer 
type, going around Zi and Z2- The correlator is then pro- 
portional to the integration from Z2 to Zi, and is written 
as 

h = M(l - ZlZlTil - ZlZ2)\l - Z2Z1Y 
x{l-Z2Z2)\zi-Z2)\zi-Z2T 

X / ' dv{l - VZifil - VZ2)\zi - V)\V - Z2f. 

J Z2 

(82) 

Here, A/j is a constant arising from the deformation of 
the contour. Note that, at this point, the expression is 
similar (in fact, proportional) to the integral represen- 
tation of a chiral four point conformal block [21] (with- 
out boundary). One may then proceed in the standard 
manner, namely, by fixing the projective SL{2, C) gauge, 
performing the integration, and then recovering the co- 
ordinate dependence. We thus have. 



h = M(i - zizini - Z2Z2TH11 - 1)]" 

r(l-a2)2 



r(2 - 2a? 



-F{2a,l - at,2 - 2at;ri), 



(83) 



where F = 2F1 is the hypergeometric function of the 
Gaussian type, and 77 is defined as 



jzi - Z2){Z2 - Zl) 

(1 - 2;i^i)(l - 2:2^2)' 



(84) 



The calculation for (II) goes similarly. As we have one 
screening operator 



Q-= 6 dvV-{v), 



(85) 



in the antiholomorphic sector, the correlator is written 
as 



/„ = u{B{ai,s)\ j> dvVi,2{zi)Vp>-i,p-2{zi)V-{v) 

X V'l,2(22)Vl,2(22)|0,0;ao)c/ 
= Ml{zi - Z2)''{Z1 - ^2)"(1 - ZiZ^r 

X(l- 0102)^(1 -^2^l)»(l-^2^2)^ 

X / ' dv{zi - vy{v - Z2f{l - ZIV)^{1 - Z2vf. 

J Z2 



(86) 
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We have again deformed Felder's integration contour 
(this time in the antiholomorphic sector) into the 
Pochhammer contour around Zi and Z2- The resulting 

integral is proportional to the one from Z2 to zi , and Afu 
is a constant. Performing the integration we have 



In = Mi(l - z,z,r{l - z2Z2T[r,{r, - 
r(l-a^_)r(3a^_-l) 



b—a 



(87) 



The two correlators /j and In with fixed boundary 
charges arc represented schematically (in the Schottky 
double picture) in Fig.l. They correspond to the two 
conformal blocks of the chiral four point function. 



D. Correlation functions on the half plane 

Boundary correlation functions obtained on the unit 
disk are mapped on the half plane by the global conformal 
transformation, 



w = -lyo 



z-l 



w 



z-l 

J+T' 



(88) 



which takes the unit circle \z\ ~ 1 on the z-plane to the 
infinite line Imw = on the w-plane, and the origin z = 
to the point w = iyo, yo G M. Under this transformation 
the holomorphic coordinate dependence on the unit disk 
is mapped on to the upper half tw-plane, and the anti- 
holomorphic dependence is mapped on to the lower half 
w-plane. Using the transformation (88), p-point correla- 
tion functions on the half plane are written using those 
on the ^;-plane, i.e. on the disk, as 



i)i(wi, wi) • • ■ (f>p{wp,Wp)) UHP 
dw 

dzi I \ dz. 



= ri(^) (^) ('^i(^i'^i)---?!'p(2;p,2;p))disk 

= n| (,,+if(°^,+i) | iMzi,z,) ■ • •<^p(.p,^p))disk, 

(89) 

where we have assumed hi = hi. The parameter tj of (84) 
is mapped under this transformation as 

^ ^ (2:1 -^2) (^2 -^l) ^ {Wl - W2){Wl - ,gQX 

(1 - ZiZi){l - Z2Z2) {Wi - Wi){W2 - W2) ' 

which is an anharmonic ratio of the four points w\, W2, 
Wl and W2- 

Now, the boundary one point function of (j)r,s on the 
upper half plane is easily found by using (70) as 



{(t>r,s{w,w)) 



_ 2?yo(l - zz) 
^("--^ " ^ (z + l)(z + l) 
-i{w-w)]-^^ 



-2h 



-2h 



(91) 



where h = hr^s = c«r,s(ar,s — 2ao) is the conformal dimen- 
sion of the operator ^r,s, and w = x+iy, iB = w* = x—iy, 
X, y G M. Note that the result is j/o-independent. The two 
point function of (/)i,2 on the disk is mapped on to the 
half plane similarly. For the conformal block (I) we have 



{4'l,2iwi,Wi)(l)i^2iw2, W2)}B{ai,i) 

2 ^ -2h 



4%^ 



{Z1 + 1){Z1 + 1){Z2 + 1){Z2 + 1) 

{Wi - Wl)(w2 - W2) 



h 



{Wi - W2){Wi - W2){W\ - W2){Wi - W2) 

r(2-2ag.) ^^~^^' 1 - ' 2 - 2«- ; ^) , 

and for (II) we have 

(<?i'l,2(M'l, W>l)^l,2(w2, W)2))B(ai,3) 



2h 



(92) 



-2h 



■Mil 



(Z1 + 1)(Z1 + 1)(Z2 + 1)(Z2 + 1) 

{Wl — W\){W2 — W2) 



2h 



{wi - W2){wi - 'W2){wi - 'W2){'Wi - W2) 

r(l-a^)r(3a^ - 1) 



'.F{al,l- a^_,2a^_;r]), 



-i-v) 



2h+2at 



(93) 



where h = hi^2 = cn,2(cn.2 ~ 2q!o)- Physical correlation 
functions are linear sums of these conformal blocks where 
the coefficients are given by Cardy's states. 



E. Boundary states and conformal blocks 

Before illustrating in specific examples, we mention 
how the above description of boundary correlation func- 
tions fits into the conventional discussion of [7,8], and see 
the validity and limitation of the Coulomb-gas approach. 

The charge neutrality conditions (67) (68) for a one- 
point function pick up the coefficient of the corresponding 
Ishibashi state from a Cardy state, since a = a^-.s is the 
only boundary charge which gives a non- vanishing term. 
This agrees with our understanding that the coefficients 
of Cardy's state are essentially the one-point coupling 
constants of bulk (closed string vertex) operators to the 
boundary (brane) [7,8]. Once one-point coupling con- 
stants are known, it is in principle possible to compute 
boundary multipoint functions since they reduce to one 
point functions after repeated use of bulk OPEs (Fig. 2a). 
In particular, we may start such a procedure from the 
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farthest point from the boundary (in the radial ordering 
sense), approaching the boundary by performing OPE 
with the farthest remaining point one by one. Due to 
naturahty of CFT, such OPEs are translated into the 
fusions of operators, 

[<t>i] X m = Ln], 
b'l] X [03] = [32], 

[ip-2] X [(Ap] = [ip-i], (94) 

which define a conformal block with no subchains 
(Fig. 2b). As the fusion of primary operators is equivalent 
to the map (55) between CBFSs restricted to BRST sub- 
spaces [14], the conformal block of Fig. 2b is represented 
by our fixed boundary-charge correlator (61) with a cor- 
responding to [ip-i]. The Ishibashi state acts as 
a filter (or a half mirror) transmitting only the Virasoro 
representation [jp-i]. This property of Ishibashi states is 
captured by boundary charges and charge neutrality. 

Due to the absence of internal channels, two point func- 
tions on the disk (or half plane) are completely deter- 
mined by a boundary state, apart from the normalisation 
of conformal blocks. In the case of p-point functions with 
p > 3, however, states in the internal channels ji<p-i 
cannot in general be determined uniquely even if the state 
in the 'boundary channel' jp-i is fixed (an example is the 
spin three point function of the Ising model with j2 = a, 
where ji can be / or e). Corresponding to this, the con- 
tours of (61) with p > 3 may be deformed in several 
different ways to give independent convergent functions 
which are expected to reproduce conformal blocks with 
different internal states. The relative coefficients of such 
conformal blocks in a boundary correlation function can- 
not be determined by the boundary state (as these con- 
formal blocks belong to the same Ishibashi state) but 
should be constrained by information of the bulk. 



(a) OPE picture (b) Conformal blocks 




FIG. 2. OPE and conformal blocks of boundary p-point 
functions. Repeating OPEs in the bulk, the boundary corre- 
lation function reduces to one point functions (a). This can 
also be seen as conformal blocks with internal channels ji<p-i 
and boundary channel jp-i (b). 



IV. ISING MODEL 

We shall illustrate the method presented in the pre- 
vious section in the example of the critical Ising model. 
Before starting actual calculations we notice from the 
discussions of the previous section that the non-unitary 
sector and the boundary states with charges outside the 
Kac table do not contribute to correlation functions. Ne- 
glecting such unnecessary terms in (51) - (53), for the 
bra-boundary states of the Ising model we have, 

(7| ~2-l(t;(B(ai,i)| + c;(B(a2,3)| 

+[/(B(a2,i)| + c/(B(ai,3)l) 
+2-3/4(^(B(ai,2)| + u{B{a2,2)\), (95) 

(e| -2-i(c/(B(ai,i)| + a(B(a2,3)| 
+c/(S(a2.i)| + [/(B(ai,3)|) 
-2-^'\u{B{ai,2)\ + u{B{a2,2)\), (96) 

(^1 ~ 2-V2(^(B(ai,i)| + c/(S(a2,3)l 

-u{B{a2,i)\-u{B{a^,i)\). (97) 

Since these are linear combinations of fixed boundary- 
charge states jj{B{ar_s)\, the correlation functions on the 
disk with physical boundary conditions /, e and a are 
given by linear combinations of fixed boundary-charge 
correlators (61). Using the global conformal transforma- 
tion explained in the last section, we shall obtain correla- 
tion functions on the half plane and compare them with 
existing results. 

A. One point functions of spin and energy operators 

Let us first consider the spin one point function. We 

may choose either a = (61.2 or cr = (f>2.2- Let, for def- 
initeness, (r, s) = (1, 2) in the holomorphic part and in 
order to use the result of the last section, (f, s) = (2, 2) 
in the antiholomorphic part. From (70) we immediately 
notice that only the state u{B{ai^2)\ contributes to the 
one point function, all other states giving vanishing cor- 
relators. For the boundary condition I, the one point 
function on the disk is 

(/>(z,z)|0) 

= (/|yi,2(z)V'2,2(^)|0,0;ao)c/ 

= 2-3/V(.B(ai,2)|Vi,2(2)t^2,2(^)lO, 0; ao)u 

= 2-3/4(1 _ zz)-'/^. (98) 

Properly normalised one point function is then, 

iMlllM = 2V4(l - ..-)-V8. (99) 
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This is mapped on to the half plane by using the confor- 

mal transformation (91), as 



(cr(usw'))/ = {<y{y))i 



,1/4 



(2?y) 



-1/8 



(100) 



where y is the distance from the boundary. Likewise, 
boundary spin correlation functions for the conditions e 
and a are obtained simply by picking up the coefficients 
of c/(i?(ai,2)| in (96) and (97), and are normalised using 
(e|0) = 1/2 and (a|0) = l/s/2. On the half plane, they 



are 



and 



{cj{w, w))~, = {a{y)), = -2V4(2y)- V8, (loi) 



{a{w,w))^ = {a{y))^=Q. (102) 

Hence, / and e are indeed the fixed (up and down) and 
a is the free boundary condition, as is stated in [6,7]. 
In our Coulomb-gas formalism the relation between one 
point functions and the coefficients of Ishibashi states is 
explained by the neutrality of charges. 

Putting, say, (r, s) = (2,1) and (f , s) = (1,3), the 
energy one point function is obtained similarly. On the 
half plane we have, 



{<y))i={e{y)h^{2y)-\ 
{e{y)), = -{2y)-\ 



(103) 
(104) 



B. Spin two point function 

Next, let us consider the spin two point function. Since 
= (l>i,2, we can use the result of Subsec.C of the pre- 
ceding section. There are only two values of boundary 
charges, — ai^i and — Q;i,3, which give non-trivial contri- 
butions to the correlator. The two corresponding states 
u{B{a\^x)\ and u{B{ax^z)\ give rise to the two conformal 
blocks Ji and In, respectively, and the correlation func- 
tion is a linear combination of these conformal blocks 
with coefiicients given by Cardy's states (95) - (97). 
Then, 

{i\a{zi,zi)(j{z2,Z2)\0) {e\a{zi,Zi)a{z2,Z2)\Q) 



m 



= h + hi, 



{a\a{zi, zi)(t{z2, Z2)\(i) 



h - hi, 



(105) 
(106) 



where the actual forms of h and hi are given by (83) and 
(87), with a = -1/8, 6 = 3/8 and = 3/4. In this case 
the hypergeometric functions reduce to, 



2^1 (- 



1 1 1 

'4'4' 2 



V2 



3 1 3 V-2(l-vT^ 
'^i^4'4'2'^^ ^ 



(107) 



(108) 



Using (107), (108) and the conformal transformation 
(88), wc find the two point functions on the half plane, 

{a{wi,wi)a{w2,W2)) i 
= {a{wi,wi)a{w2,W2))e 

_ \ {Wi - Wi){W2 - W2) 1^^^ 

■v/2 \ (Wl - 'W2){W1 - W2){wi - W2){WI - W2) j 

X (^i^ +1 + Mi/T^ - 1) . (109) 

{a{wi,Wi)(T{w2,W2))a 

_ ( jwi - Wl)(w2 - W2) 1 

■v/2 \ (Wl - W2){wi - iU2){'Wi - iU2){'Wl - W2) j 

X (mi^/^/T^^+I - Afu\/V^^V^^ , (110) 

where J\f = r(l/4)Vr(l/2). Studying the behaviours 
away from the boundary (we accept the convention that 
two point functions of bulk operators are normalised 
as {(f)i{wi)(f)j{w2)) = Sij{wi — ^2)"^''') and compar- 
ing the leading terms with the OPE coefficients of [8], 
the normalisation of the conformal blocks are fixed as 

M = Mi = r(i/2)/r(i/4)2. 

This result was obtained long time ago [13], by solv- 
ing a differential equation to find the two conformal 
blocks h, hi, and then fixing the coefficients by con- 
sidering asymptotic behaviours of the correlation func- 
tion. The (relative) coefficients of the conformal blocks 
are now attributed to the coefficients of Cardy's states in 
our Coulomb-gas approach, although we have used the 
asymptotic behaviours to fix the normalisation of each 
conformal block. 



C. Energy two point function 

Finally we derive the energy two point function in our 
formalism. As e = (^2,1, the calculation is parallel to the 
case of the spin two point function. Prom the charge neu- 
trality condition we find that non-vanishing correlators 
arise from the states u{B{ai^i)\ and u{B{a3^i)\. How- 
ever, as none of the boundary states (95) - (97) contain 
u{B{a^_i)\, only u{B{ai_i)\ gives non-trivial contribu- 
tion to the correlation function. Hence, the energy two 
point function does not depend on boundary conditions. 
After a simple calculation we find, on the half plane, 



{e{wi,w-i)e{w2,W2)) i^^^g. 



{wi — Wi){W2 — W2) 



{Wi - W2)[wx - W2)(W\ - W2){W\ - W2) 

?(-2/3)' ^(-^' - -2/3; n). (Ill) 

The normalisation constant M is determined as A/" = 
r(— 2/3)/r(— 1/3)^, by considering the off-boundary be- 
haviour. The hypergeometric function turns out to be an 
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algebraic function F(— 2, —1/3, —2/3; ry) = 1 — ry + r/^, and 
using the coordinates Wi = Xi + iyi, Wi = w* = Xi — iyi 
{xi, Hi e M), the correlation function (111) is written as 

(e(a;i,yi)e(a;2,2/2))/,e,a 

- x^y + {yi - y2y][{x^ - x^Y + (yi + ys)^] 



This agrees with the result in [13]. 

V. SUMMARY 

In this paper we have described a novel method to cal- 
culate correlation functions of two dimensional CFT on 
the disk and on the half plane. We have used the free-field 
construction of boundary states developed in [18], and 
derived boundary correlation functions for the bound- 
ary states classified by Cardy's method. The key feature 
of our formalism is the neutrality of bulk and bound- 
ary charges, which associates the coefficients in Cardy's 
boundary states directly with the linear combinations of 
conformal blocks. Thus we could unify the two important 
pieces of boundary CFT, namely, boundary correlation 
functions [13] and consistent boundary states of Cardy 
[6], by using the Coulomb-gas picture. We have checked 
the formalism in the Ising model, and shown that our 
method reproduces the known results. 

Cardy's classification of boundary states has been gen- 
eralised by Lewellen [8] and Pradisi, Sagnotti, Stanev 
[9] beyond the diagonal models, and the Coulomb-gas 
technique is also known to be applicable to more general 
CFTs, such as WZNW models [29] and CFTs with W- 
algebra [30-32]. We therefore expect that the method 
discussed in this paper may be applied to such CFTs rel- 
atively easily. In particular, from a string theory point of 
view, application to WZNW theories seems to be quite 
fruitful since it would provide an alternative method to 
find correlation functions with D-branes on a group mani- 
fold. We shall discuss such issues in separate publications 
[33]. 
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